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Ultrarelativistic electron-positron plasmas can be produced in high-intensity laser fields and play 
a role in various astrophysical situations. Their properties can be calculated using QED at finite 
temperature. Here we will use perturbative QED at finite temperature for calculating various im- 
portant properties, such as the equation of state, dispersion relations of collective plasma modes 
of photons and electrons, Debye screening, damping rates, mean free paths, collision times, trans- 
port coefficients, and particle production rates, of ultrarelativistic electron-positron plasmas. In 
particular, we will focus on electron-positron plasmas produced with ultra-strong lasers. 
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lO ■ Ultrarelativistic plasmas, in which the thermal energy of the particles is much larger than their rest mass energy, 
were discussed first in the fifties in the context of astrophysics. They occur in the early Universe as well as stellar 
and galactic high-energy processes (Zel'dovich and Novikov, 1971; Zel'dovich and Novikov, 1983; Raffelt, 1996). For 
example, ultrarelativistic electron-positron plasmas (EEP) can be created either by high temperatures in supernovae 
[ explosions (see e.g. Hardy and Thoma (2001)) or strong magnetic fields in so-called magnetars (Beskin et ai, 1993). 
■ The theoretical description of ultrarelativistic plasmas is based on transport theory (Silin, 1960) or thermal field 
' theory (Tsytovich, 1961). The latter is an extension of quantum field theory, such as QED or QCD, in vacuum to 
^ , finite temperature (or chemical potential) describing high-energy particle interactions in matter. For this purpose, 
' two different approaches were developed, the imaginary time formalism (Matsubara, 1955; Kapusta, 1989) and the 
ijjj , real time formalism (Landsmann and van Weert, 1987). 

' In the eighties a new type of ultrarelativistic plasmas became of interest, the quark-gluon plasma (QGP) (Miiller, 
1985). It corresponds to a new state of matter of deconfined quarks and gluons at extremely high temperatures 
above ksT = 150 MeV. The early Universe should have been in this state for the first few microseconds after the Big 
Bang. Furthermore, the QGP is expected to be created in relativistic heavy-ion collisions, in particular within the 
accelerator experiment RHIC (Gyulassy and McLerran, 2005). There a hot and dense fireball ("little bang") of the 
size of an atomic nucleus is produced, which could be in the QGP phase for less than 10~^^ s. Therefore the QGP 
in those experiments cannot be observed directly but its discovery relies on the comparison of theoretically predicted 
signatures with experimental data ("circumstantial evidence"). For understanding the properties of the QGP and 
for calculating signatures for its formation, QCD at finite temperature has been applied. Besides a non-perturbative 
method based on a numerical solution of the QCD equations (lattice QCD), perturbation theory at finite temperature 
has been used widely (Thoma, 1995a; Le Bellac, 1996). Lattice QCD allows only the calculation of static quantities, 
and hence not of most of the proposed signatures following from particle production. Perturbative QCD, on the other 
hand, has been used for computing static as well as dynamic quantities. However, its predictions are limited by the 
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fact that the QGP is a strongly-coupled plasma at temperatures achievable in the experiments. (At extremely high 
temperatures far above the transition temperature perturbative QCD should be reliable due to asymptotic freedom 
of QCD.) Nonetheless, a wide variety of properties and quantities of the QGP have been considered in this way. It 
started with the calculation of self-energies and dispersion relations in the high-temperature limit (Klimov, 1982; 
Weldon, 1982a; Weldon, 1982b). In the high- temperature limit the results differ from the QED results only by trivial 
factors containing the QCD degrees of freedom (color, flavor). In addition, it can be shown that the gluon self-energy 
or polarization tensor, which is directly related to the dielectric function (see e.g. Elze and Heinz (1989)), can be 
derived also within the transport (Vlasov) approach (Silin, 1960), since the high-temperature limit corresponds to the 
classical limit. Other quantities, however, such as damping or production rates and transport coefficients, require a 
quantum field theoretic treatment. In particular, a resummation of certain diagrams, the hard thermal loops (HTL), 
is needed for a consistent description, i.e. to obtain gauge invariant, infrared finite (or improved) results complete to 
leading order in the coupling constant (Braaten and Pisarski, 1990). Details about quantum field theoretic methods 
and their application to the physics of the QGP can be found in the review article by Thoma (1995a). 

Recently a third possibility to study ultrarelativistic plasmas has been suggested: in extremely strong laser fields 
the creation of an ultrarelativistic EPP with temperatures around 10 MeV could be realized soon (Liang et al., 
1998). For example two opposite, circular polarized laser pulses with a duration of 330 fs and an intensity of 7 x 10^^ 
W/cm^ could hit a thin gold foil. In this way the target electrons can be heated up to about 10 MeV producing an 
ultrarelativistic EPP by pair creation. The positron density could be about 5 x 10^^ cm~^ (Shen and Meyer-ter-Vehn, 
2001). This will allow for the first time to study an ultrarelativistic EPP in the laboratory. Therefore predictions of 
the physical properties of such a system are needed. 

The interaction of relativistic electrons and positrons is described by QED. In the presence of a thermal plasma 
background (heat bath) QED at finite temperature has to be considered. As discussed above, perturbative methods 
based on the imaginary (Matsubara) or real time formalism have been developed and applied to the physics of a QGP. 
Here wc will transfer the methods and results for properties of a QGP to the case of an EPP, where a perturbative 
treatment is more reliable than in the case of a strongly coupled QGP. Some applications of thermal field theory 
to astrophysical plasmas have been discussed before (see e.g. Thoma (2002) and references therein and Altherr and 
Kraemmer (1992)). 

An important difference to non-relativistic ion-electron plasmas (Lifshitz and Pitaevskii, 1981) are the relevant 
scales. In the non-relativistic case case they are given by the masses of the plasma particles and the temperature T. 
For example, the electron plasma frequency reads 



where pe is the electron number density, the temperature of the electron component, and rrig the electron mass. 

In an ultrarelativistic plasma with T m the masses can be neglected and the important scales are the temperature 
T, called the hard scale, and the soft scale eT, which determines the collective phenomena as we will see below. Here 
we use natural units, i.e. ft = c = fc^ = 1, as usual in quantum field theory. In these units e = 0.3 corresponding to 
a fine structure constant a = e^/(47r) = 1/137. For converting to Sl-units we use 1 MeV = 1.60 x 10"-^^ J = 5.08 x 
10^2 m-i = 1.52 X 10^1 s-^ 

In the next section we will discuss the equation of state of an equilibrated ultrarelativistic EPP. Then collective 
phenomena will be considered. Afterwards transport properties and particle production will be discussed. Finally we 
will describe properties of an EPP which is not in chemical equilibrium as in the case of laser induced plasmas. We 
will not consider here the formation process and equilibration of an EPP. 

Many results presented here can be found in the literature for the case of a QGP (Thoma, 1995a) differing only 
by numerical factors, e.g. number of degrees of freedom. The purpose of this colloquium is to summarize these 
results and to extend them to laser induced EPPs as a reference for future experiments. As an example we consider 
a temperature of 10 MeV as it can be typically realized in laser produced and supernovae EPPs. 

Laser produced QED plasmas have also been discussed recently in two review articles (Mourou et ai, 2006; Marklund 
and Shukla, 2006) with emphasis on the production mechanism and non-linear effects. Here, however, we want to focus 
on the properties of an equilibrated EPP as they can be calculated from perturbative QED. Such an EPP in thermal 
and chemical equlibrium might be the outcome of future laser experiments if the intensity can be increased further on. 
As an example, we have chosen the predictions of Shen and Meyer-ter-Vehn (2001) based on a numerical simulation 




(1) 



and the Debye screening length due to the electrons in the plasma 
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(PIC) and cross sections for electron-positron productions, in which two opposite laser beams are focussed on a thin 
gold foil leading to a chemically non-equiliibrated plasma (see section VI). However, this interesting proposal still 
waits for its experimental confirmation. Other production mechanisms, based on the Schwinger pair production effect 
in strong fields (Schwinger, 1951), have shown that pair production can be efficient already far below the Schwinger 
critical field strength, requiring laser intensities of 5 x 10^^ W/cm^, in the case of time dependent and inhomogeneous 
fields, e.g. two oppositely directed pulsed laser beams in vacuum (see e.g. Avetissian et al, 2002; Narozhny et al., 
2004; Di Piazza, 2004; Dunne and Schubert, 2005; Gies and KlingmuUer, 2005; Blaschke et al. 2006a, Blaschke et al., 
2006b, Schiitzhold et al, 2008). Furthermore the pair production in an X-ray free electron laser (XFEL) has been 
discussed (Ringwald, 2001; Alkofer et al., 2001; Roberts et al., 2002). QED plasmas in strong magnetic fields have 
been considered by Danielsson and Grasso (1995) and more general in strong electromagnetic fields by Morozov et al. 
(2002). QED plasmas have also been studied in great detail in the book by Melrose (2008). 



II. EQUATION OF STATE 

In this section we will consider the equation of state of an EPP under the following assumptions: 

1. ultrarelativistic EPP, i.e. T ^ m, 

2. thermal and chemical equilibrium, 

3. equal electron and positron density, i.e. vanishing chemical potential, 

4. ideal gas, i.e. no interactions in plasma, 

5. infinitely extended, homogeneous, and isotropic EPP. 

We will relax some of these assumptions in the following sections. According to these assumptions the distribution 
function of the electrons and positrons is given by the Fermi-Dirac distribution 

and of the photons by the Bose-Einstein distribution 

npip) = ^p/T^^ , (4) 

where the momentum p is identical to the energy E of the particles in the ultrarelativistic case. It should be noted 
that the photons are in equilibrium with electrons and positrons under the above assumptions, i.e. the system is 
actually an electron-positron-photon gas. 

The particle and energy density can be calculated by integrating over the distribution functions. In the case of a 
QGP the results can be found in Miillcr (1985). The particle number density of the electrons and positrons follows 
from integrating over the Fermi-Dirac distribution as 

= / (S^"^^^^ = ^ ^^^^ = 

where gp = 4 is the number of degrees of freedom corresponding to the electrons and positrons in the two spin states. 
Assuming a temperature of T = 10 MeV, we find pi" = 370 MeV^ = 4.9 x 10^° m-^. 

The photon density follows accordingly by integrating over the Bosc-Einstcin distribution with gs = 2 degrees of 
freedom corresponding to the two polarization states as = ({3) = 0.24 T^. 

The energy density of the electron-positron-photon gas is obtained from 

=''j (0^"^(^^ + / (0^^"^(^^ = ^ = ^-'^ 

where the photons contribute 36 % to the energy density. Here the Boltzmann law, e*^' ~ T^, holds also for the 
fermions because we neglected their masses. 

For T = 10 MeV we find e*^* = 3.8 x 10^^ J m^^. In a volume of 10^^ m^ (corresponding to the size of a neutron star) 
the total thermal energy of the EPP is 3.8 x 10*^ J, which corresponds to about 10% of the entire energy (without 
neutrinos) released in a supernova type II explosion. In a volume of 1 iivo? there is still an energy of 3.8 x 10^^ J 
contained. 

The Coulomb coupling parameter of the EPP, which is a measure for the non-ideal behavior of a plasma (Ichimaru, 

1982), is given by A = e^/{dT), where d ~ Pe"^"^/^ = 2.7 x 10"" m is the interparticle distance. For T = 10 MeV we 
find A = 5.3 X 10""^ which shows that the EPP is a weakly coupled plasma in contrast to the QGP where A = 0(1) 
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(Thoma, 2005) . Therefore the ideal gas results for the equation of state derived above are a good approximation. After 
all, interactions in the EPP play an important role, for example, for the collective behavior of the plasma as discussed 
in the next section and for equilibration of the plasma. Obviously, the interaction can be treated by perturbation 
theory. 



III. COLLECTIVE PHENOMENA 



Interactions between particles can be separated into two classes: individual collisions between the particles and 
long-range interactions of particles with the medium (Lifshitz and Pitaevskii, 1981). The latter ones lead to collective 
effects, which are characteristic for plasmas. The crucial quantity from which the collective phenomena are derived is 
the dielectric tensor relating the macroscopic electric field Di in the medium to the external field Ei {i — x, y, z), i.e. 
in momentum space 

A(c^,k) =^ey(c^,k)i;,(c^,k). (7) 
j 

In the case of an isotropic medium it depends only on fc = |k| and has two independent components 

e.,{u,k) = eTico,k) i^5,,-^^+e,{co,k)^ (8) 

The dielectric tensor is closely related to the polarization tensor or photon self-energy by (see e.g. Elze and Heinz 
(1989)) 



nL(cj,fc) 

eL(uj,k) = 1 — 

nT{LJ,k) 



eriuj.k) = 1 - 



(9) 



where 11/, and IIt are the longitudinal and transverse components of the polarization tensor, respectively. 

The lowest order diagram for the polarization tensor is shown in Fig.l. Assuming the external momentum to be 
soft, i.e. uj and k to be much smaller than T and the internal loop momenta to be hard (HTL approximation) an 
analytic result can be found using the real or imaginary time formalism (Klimov, 1982; Weldon, 1982a) 



ni(c^,fc) 
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(10) 



where niph = eT/3 is called the effective photon mass. For T = 10 MeV we get ruph — 1 MeV. 

The hard thermal loop approximation, which is identical to the high-temperature limit, produces consistently the 
lowest order result in finite temperature QED (Braaten and Pisarski, 1990). The dielectric functions following from 
(O and (fTU]) can also be derived from the classical Vlasov equation together with the Maxwell equations (Silin, 1960), 
since the high-temperature limit corresponds to the classical limit. 
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FIG. 2: Photon dispersion relation. 

The dispersion relations of collective plasma modes, i.e. propagation of electromagnetic waves in the plasma, can 
be found by using the Maxwell equation, leading to (see e.g. Carrington et al. (2004)) 

eL{uj,k) = 0, 

erito^k) = —. (11) 

Combining ([9]), (jlOp . and (fTTj) gives the dispersion relations ujL,T{k) of the transverse as well as longitudinal plasma 
waves as shown in Fig. 2. The longitudinal branch, which does not exist in vacuum, is called plasmon as in the case of 
non-relativistic plasmas. The transverse branch does not play a role in non-relativistic plasmas but is equally important 
as the longitudinal one in relativistic plasmas. Both branches start at the plasma frequency LOpi = LOL.rik = 0) = mph 
(Kajantie and Kapusta, 1985). Consequently the collective plasma waves have soft momenta of the order eT. At 
high momenta k ^ mpii the transverse mode approaches the free dispersion ujt = k, corresponding to a real photon 
in vacuum, whereas the longitudinal mode disappears, i.e. its spectral strength is exponentially suppressed (Pisarski, 
1989). For T = 10 MeV we find iVpi = 1.5 x lO^^ s~^. 

Electromagnetic plasma waves in an ultrarelativistic EPF have also been discussed by Medvedev (1999). For the 
plasma frequency the same dependence on the temperature and coupling constant was found (see eq.(3) in Medvedev 
(1999)) if one uses the equilibrium number density for the photons proportional to (see above) there. However, 
the numerical pre-factor there is wrong and only one branch of plasma waves was discussed. 

Another important quantity which can be derived from the polarization or dielectric tensor is the Debye screening 
length, entering the Yukawa potential of a heavy, non-relativistic test charge in the EPP. The Debye screening length 
is given by the static limit of the longitudinal component of the polarization tensor 1/11^(0; = 0) (Kajantie and 
Kapusta, 1985), leading to Ad = l/{Vimph), which is 1.1 x lO^^^ ni for T = 10 MeV. 

Finally from PU)) we see that the polarization tensor and the dielectric function become imaginary for uj'^ < k^, 
i.e. below the light cone uj = k, corresponding to Landau damping (Pisarski, 1988). We also observe that the plasma 
waves calculated at lowest order perturbation theory are undamped since they are located at w > fc. This changes, 
however, at higher orders where the dispersion relation can intersect the light cone (Carrington et ai, 2004). It is 
interesting to note that all these results are apart from some numerical color and flavor factors also valid for collective 
gluon modes in the QGP in the high-temperature limit (Thoma, 1995a). 

A complete new phenomenon that does not appear in non-relativistic plasmas is the existence of fermionic plasma 
waves because all fermion masses are much too large in the non-relativistic case. Their dispersion relations follow 
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from the pole of the electron propagator containing the electron self-energy. Using again the hard thermal loop 
approximation the electron self-energy reads (P = {pq,p), p = |p|) (Klimov, 1982; Weldon, 1982b) 

E(P) = -a{po,p) P^j^ - bipo,p) 70 (12) 

with 



a{Po,p) = ^ [ir(P''7^S)-poir-(7oS)] 



b{po,p) = ^ [P^ tr(7o S) - po ^r(P^7M S)] , (13) 

where the traces over the 7 matrices are given by 

triP'^j^i:) = Ami, 

fr(7oS) = 2m%- In^^^ (14) 

P Pa-P 

with the effective electron mass mp = eT/y/S (Klimov, 1982; Weldon, 1982b), which is 1.1 MeV at T = 10 MeV. 
The full electron propagator in the helicity representation is given by (Braaten et al., 1990) 

= 2D^) - ^ • ^) + ^D^P) i^o+P-l), (15) 

where 

D±{P) = -po±P+^ [±tr{P'^l^ S) - {±Po - P) Hjo S)] . (16) 

Again these results agree with the ones for quarks in the QGP apart from numerical pre-factors. 

The dispersion relations following from the pole of this propagator are shown in Fig. 3. Again two branches show 
up, one with a positive ratio of the helicity to chirality following from = 0, the other one with a negative 
ratio (ui-) following from £>_ = 0, called plasmino (Braaten et al., 1990). The plasmino branch w_, which does not 
exist in vacuum, shows an interesting behavior, namely a minimum at fc = 0.41mF. The same behavior of the quark 
dispersion in the QGP has been foimd. In this case sharp peaks in the dilepton production rate can appear due to the 
minimum in the plasmino branch, which leads to van Hove singularities (Braaten et al., 1990; Peshier and Thoma, 
2000). Whether something similar could be observed in the EPP, e.g. in the electron spectrum, is a very interesting 
question which should be investigated in detail. It could open the exciting possibility to observe a new collective 
plasma wave, the plasmino, experimentally in a laser induced EPP. 

The collective qiiantities derived here perturbatively from the dielectric functions are linear phenomena. Nonlinear 
collective effects in electron-positron plasmas, have been considered in a number of publications, e.g. soliton formation 
(Lontano et al, 2001), nonlinear Alfven waves (Zhao et al., 1994), nonlinear photon interactions (Tajima and Taniuti, 
1990), plasma enhanced photon splitting (Brodin et al., 2007), and nonlinear self-modulation of radio pulses (Chian 
and Kennel, 1983). 



IV. TRANSPORT PROPERTIES 



Now we want to consider the interaction and properties of particles in the plasma with hard momenta, i.e. of 
the order of T or larger. In particular we are interested in damping and transport rates, mean free paths, collision 
times, energy losses of these particles and other transport properties such as the shear viscosity of the EPP. All 
these quantities have been calculated perturbatively in high- temperature QCD (Thoma, 1995a). As in the case of the 
collective phenomena, the results can be almost directly taken from the QCD calculations. 

It was shown by Braaten and Pisarski (1990) that a consistent treatment of gauge theories such as QED at finite 
temperature, i.e. for obtaining results that are gauge independent, infrared finite, and complete to leading order, 
require the use of an effective perturbation theory using resummed Green functions based on the HTL approximation 
(HTL resummation technique) . The HTL method relies on the assumption that one can distinguish between soft and 
hard momenta, i.e. eT <C T. This is approximately fulfilled for QED, where e = 0.3, but not in QCD, where the 
corresponding coupling constant g > 1. 

The damping rate of an electron or positron in the EPP is defined as the imaginary part of the dispersion relation 
^l,t{p)- To lowest order it follows from the elastic scattering diagram of Fig. 4. In the case of a hard electron or 
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FIG. 3: Electron dispersion relation. 
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FIG. 4: Lowest order diagram for electron-electron scattering. 



positron with momenta of the order of T or higher it exhibits a quadratic infrared (IR) divergence which can be 
reduced to a logarithmic one using a HTL resummed photon propagator. This logarithmic singularity is expected to 
be cut-off by higher order contributions leading to (Thoma, 1995a) 



7e 



e^T, 1 
- — In - 

47r e 



(17) 



within logarithmic accuracy, i.e. the constant under the logarithm is not determined. For T = 10 MeV we obtain 
7e — 86 keV, which is much smaller than LOpi = 1 MeV, showing that the EPP is not overdamped. 

Physically more important are the transport rates F which are related to the mean free path and collision time 
of electrons and positrons in the EPP. They differ from the damping rate in cutting off the long range interactions 
with small scattering angles by a factor (1 — cos6') under the integral defining the rate (Lifshitz and Pitaevskii, 
1981). This leads to an improvement of the IR behavior (logarithmic instead of quadratic singularity in perturbation 
theory) and a finite result using the HTL method. Logarithmic divergent quantities can be treated consistently by 
splitting them into a soft part and a hard part, where the soft part is calculated using the HTL resummation technique 
(Braaten and Yuan, 1991). For the transport rate we find to logarithmic accuracy by extending the QCD results in 




FIG. 5: Diagrams defining the collisional energy loss. 



Thoma (1994) to QED 



re = ^ln-, (18) 



where the Mandelstam variable s = {P + K)'^ is the square of the sum of the four momenta of the incoming particles 
in the scattering diagram of Fig. 4. For deriving the above result we replaced in (15) of Thoma (1994) Ca9^ by and 
•nig by mph- For going beyond the logarithmic approximation, which is valid within about a factor of two, one has to 
calculate explicitly the hard contribution from the QED scattering diagrams as it was done by Thoma (1994) for the 
QCD case. For thermal particles we replace s by its thermal average (s) = 2{p)e{k)e — 19. ST^ (Thoma, 1994), where 
(p)g = = eli/pfJ = 3.11T. Assuming again T = 10 MeV, we get = 0.54 keV. 

The mean free path A™^^ and collision time Te of the plasma particles (electrons and positrons) are given by the 
inverse of the transport rate l/Fe, leading to A™^^ = 0.37 nm and Te = 1.2 x 10"^^ s at T = 10 MeV. 

In a non-relativistic plasma the shear viscosity can be estimated from elementary kinetic theory as (Reif, 1965) 



1 ^p,(p,)Ar^p (19) 



where the sum is performed over the various components of the system. In an relativistic plasma the coefhcient 1/3 
should be replaced by 4/15 (de Groot et ai, 1980). Using the mean free path following from (fTS)) . the density of (O, 
and the thermal momentum (p)e = 3.1 IT, the shear viscosity is given by (within logarithmic accuracy) 

At T = 10 MeV the shear viscosity coefficient is rje = 7.9 x 10^" Pa s. 

Another quantity of interest in a plasma is its stopping power or the energy loss of an energetic particle in the 
plasma. There are two contributions, namely the energy loss by collisions and the radiative one by bremsstrahlung. 
In a relativistic plasma the latter one becomes important. The relevance of these contributions in the QGP has been 
discussed controversially (Mustafa and Thoma, 2005). The collisional energy loss is given by the mean energy transfer 
divided by the mean free path leading to (Braaten and Thoma, 1991) 



dE _ 1 
dx V 



d-fuj (21) 



where v is the particle velocity, 7 the damping or interaction rate proportional to the plasma density and the collision 
cross section, and w the energy transfer from the energetic particle to plasma particle in the collision. Using for the 
collision cross section the lowest order diagrams in Fig. 5, the collisional energy of a muon with mass M in an EPF 
has been calculated by Braaten and Thoma (1991) applying the HTL resummation technique 

dE e^T^ (l l-v^ \+v\ (E \ ^, \ 
dx 247r \v 2v^ 1-vJ \ M e J 

where A{v) is a slowly varying function of the muon velocity v between 1.3 and 1.5. 

The collisional energy loss of an electron with energy E T is approximately given by (Braaten and Thoma, 1991) 

dE e^T^ , 15.3E 



dx 487r e^T 



111 ^2^- (23) 



In Braaten and Thoma (1991) the factor 15.3 was replaced by 7.6 due to an numerical error which was corrected in 
Thoma (1995a). This leads to an energy loss of 200 MeV/nm for an electron or positron with an energy of E = 100 
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FIG. 6: Diagrams defining the photon damping rate. 

MeV at r = 10 MeV, showing that such an electron is stopped (thermaUzed) within a fraction of a nanometer. (Note 
that the plasma density, proportional to T^, on which the energy loss depends via the mean free path is hidden in 
the temperature dependence of the energy loss.) Recently it has been shown that this calculation can be improved by 
taking into account additional diagrams changing this result slightly (Peigne and Peshier, 2008). So far no calculations 
of the radiative energy loss in an EPP have been performed to our knowledge. 

The damping rate of a photon in an EPP follows from the diagram in Fig. 6, where a HTL resummed electron 
propagator has to be used in case of soft momenta of the exchanged electron (positron). In contrast to the electron 
damping rate, the photon rate is infrared finite using the HTL method due to the presence of an electron propagator 
in Fig. 6 instead of the photon propagator in Fig. 4. Hence there is no need to cut off the long range interaction 
introducing a transport cross section. The result for a photon with energy E = p reads (Thoma, 1995b) 

e'^T^ , 3.88£: 

The mean free path and the collision time of photons in an EPP are given by l/Tph- For a thermal photon with the 
mean momentum {p)ph — ^p'h/Pp'i ~ 2.75T at T = 10 MeV the mean free path X™/^ = 0.28 nm and the collision 
time Tph = 9.4 x 10~^^ s follow. Actually the damping rate given in (|24|) is a lower limit as higher order effects will 
enlarge it. As a matter of fact, the photon production rate in a QGP, which is the inverse process of the damping 
rate (Thoma, 1995b), was shown to be about a factor of 2 larger taking bremsstrahlung into account (Arnold et ai, 
2001). 

For the viscosity of the photon component using the above mean free path, the photon density (see above), and the 
mean photon energy {p)ph — 2.75T we find 

_ 48.7r^ 
"^"^ ~ e4ln(3.27/e) 

corresponding to 3.5 x 10^° Pa s at T = 10 MeV. Hence the viscosity of the EPP ?y = rje+Tjph has similar contributions 
from the electrons and photons. 

A more advanced calculation of the total viscosity of the EPP based on the Kubo formula yields within logarithmic 
accuracy (Arnold et ai, 2000) 

188T3 ^ ^ 

This result is about a factor of 1.5 larger than the one presented here based on the elementary kinetic theory, which 
is typically valid within a factor of 2 (Reif, 1965). 

V. PARTICLE PRODUCTION RATES 

At high temperatures above 10 MeV also other particle species will be produced, e.g. muons with a mass of nip = 106 
MeV. Their rate follows to lowest order from the diagram in Fig. 7 (Born term). We assume that nie ^ T ^ irip 
holds. The first inequality implies that the electron mass can be put to zero and the latter inequality implies that 
muons are not equilibrated. Then the muon production rate to lowest order (e^e^ — > 'y* — > ^^/i+) is given by 

dN _ (^^2ml\(^ ^^iX'^ T 1 ,^ 1 + exp[-(i? + p)/(2r)] ^^7) 



d^xd^p 247r4 I AP \ AP / p exp(£:/r) - 1 1 + exp[- {E - p) / {2T)y 
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FIG. 7: Lowest order contribution of the muon production. 



where = E'^ — is the invariant mass of the virtual photon 7*, E its energy and p = |p| its momentum. This 
formula was derived by combining the cross section for the process in Berends et al. (1973) with the production 
rate in Cleymans et al. (1987) where the quark-antiquark annihilation process qq —^ fJ-^ fJ-^ was considered. Here we 
assumed the chemical potential fi — 0. The only difference to the process considered here comes from the fractional 
charge of the quarks, i.e. we use here ~ 1' ^^'^ from the number of colors in the distribution functions in 

(6a,b) of Cleymans et al. (1987), i.e. dividing the QGP rate by a factor of 9. In addition a factor 3/2 appears in 
the electron-positron annihilation cross section - see (6) of Berends et al. (1973) after integrating over the angle - 
compared to the quark-antiquark annihilation cross section in (5) of Cleymans et al. (1987), using the relative velocity 
of the ultrarelativistic particles Vqq = 1. Because of = E^ — p^ > 4:m^ the rate is suppressed exponentially for 
temperatures below 2m^. 

In order to obtain the spectrum from this formula one has to integrate over the space-time volume, taking into 
account the space-time evolution by using, for example, a hydrodynamical model. The total muon yield then follows 
from integrating the spectrum over the energy and momentum of the virtual photon. 

At temperatures above 10 MeV also hadron production becomes important, in particular pion production 
(Kuznetsova et al., 2008). 



VI. NON-EQUILIBRIUM AND FINITE CHEMICAL POTENTIAL 



EPPs produced in strong laser fields are probably not in complete equilibrium. For example, it has been predicted 
by Shen and Meyer-ter-Vehn (2001) that an positron density of about 5 x 10^* m"'^ at a temperature of 10 MeV 
can be reached. This density deviates from the equilibrium density ([5]) by 12 orders of magnitude. In the following 
we will therefore assume that the EPP produced by lasers is in thermal but not in chemical equilibrium. Then we 
can replace the distribution functions for the electrons by Fermi-Dirac distributions multiplied by a fugacity factor 
A describing the deviation from chemical equilibrium, /f(p) = ^n-p^p). This assumption has been used for example 
for describing the chemical equilibration of the QGP in ultrarelativistic heavy-ion collisions (Biro et al., 1993). The 
fugacity is given by the ratio of the experimental to equilibrium particle density, since the experimental density follows 
from integrating over the non-equilibrium distribution, i.e. 



d^p 
(2nf 



Oexp = gP / /r, .^ ^npip) = Xpeq ^ A = 10 (28) 



Using the real time formalism, QED perturbation theory and the HTL method can also be extended to non- 
equilibrium situations like the one discussed above (Carrington et al., 1999). For example, the effective photon mass 
is given now by 

4g2 /-oo 

= ^ dPPfpip) (29) 

For T = 10 MeV we then find for the non-equilibrium photon mass Wp^"^' = \/Xmph = 1 eV and the plasma 
frequency cj^""'^^ = 1.5 x 10^^ Hz. The Debye screening length in such an EPP is Ad ~ 0.1 /im. In order to speak of 
a plasma with collective behavior its dimension L should be much larger than Xjj, i.e. at least of the order of 1 ^m. 

Furthermore, an anisotropic EPP can also be described by quantum field theoretic methods (Mrowczynski and 
Thoma, 2000). In this case instabilities can occur (Mrowczynski and Thoma, 2007). 
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Finally a possible difference between the positron density and the electron density can be treated by introducing a 
finite chemical potential fi, i.e. using the distribution 

for the electrons (negative sign) and positrons (positive sign). Such a difference comes from the fact that the laser 
produced EPP is embeded in a hot electron and cold ion background of the target. Therefore there will be an excess 
of electrons over positrons in the hot, relativistic EPP. The methods described above, such as the HTL resummation, 
can be generalized easily to this case (Vija and Thoma, 1995). For example, the energy density is given by 

....li^r' + lrv + ji,/ (31) 



or the effective photon energy by 



= + (32) 



The muon production rate for finite chemical potential can be found by replacing the last three factors in (|27p by (13) 
of Cleymans et al. (1987). 



VII. CONCLUSIONS 



The properties of ultrarelativistic EPPs can be treated by QED perturbation theory at finite temperature similarly 
as the QGP by QCD using the imaginary or real time formalism. A consistent treatment requires the use of the HTL 
resummation technique for many quantities under consideration. Many results can be more or less directly taken 
from the corresponding QGP calculation. The examples discussed here are the equation of state, the dielectric tensor, 
the dispersion relations of photons and electrons, the Debye screening length, damping and transport rates, mean 
free paths, collision times, shear viscosity, stopping power, and the muon production rate. In particular collective 
phenomena, such as plasma waves, are of interest. In an EPP completely new phenomena appear, namely collective 
fermion modes associated with the possibility to observe van Hove singularities. 

EPPs produced by strong lasers open the unique possibility to investigate these plasmas and their properties which 
also exist in astrophysical systems, e.g. supernova explosions. However, current predictions indicate that laser induced 
EPPs are not in chemical equilibrium. Their properties therefore require an extension of perturbative QED and the 
HTL resummation technique to non-equilibrium situations, which is possible within the real time formalism. As a 
first application we discussed the Debye screening length in such an EPP, demonstrating that the EPP should have 
a spatial extension of at least one micron to exhibit a characteristic plasma behavior. 

Here we have considered a number of relevant properties of an ultrarelativistic, equilibrated EPP as a reference for 
future laboratory experiments and astrophysical EPPs and gave numbers for a plasma temperature of 10 MeV. The 
results are summarized in the table below. 

New results, which have not been published so far, are the electron transport rate (fT8|) . the electron viscosity (|22p . 
and the photon viscosity ([?f|) . Also new is the application of the results, derived here or compiled from the literature, 
to relativistic laser plasmas as presented in the third row of the table and in the discussion of non-equilibrium effects. 

Further quantities, e.g. the radiative energy loss, or higher order corrections, e.g. beyond the leading logarithm 
approximation, to the quantities considered above could be evaluated in a similar way as discussed above. Also an 
extension of the properties discussed in sections II - V to the case of a chemical non-equilibrated EPP would be of 
interest. Finally we did not consider here the formation of an EPP in a strong laser field, e.g. by the Schwinger 
mechanism, and its equilibration. 

Acknovifledgment: I would like to thank D. Habs, J. Rafelski, I. Kouznetsova, and G. Moore for helpful discussions 
and hints. 
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Table: Properties of an EPP 



Quantity 


Formula 


Value at T = 10 MeV 


Electron- Positron Density 


pl" = 3/7r^ C(3) 

Pll = 2/^2 C(3) 


4.9 X 10*" m"^ 


Photon Density 


3.2 X 10"*" m"^ 


Electron-Positron Energy Density 


2.4 X 10^^ J m-3 


Photon Energy Density 


t^" = llTrVeOT* 


1.4 X 10^^ J m-3 


Total Energy Density 


3.8 X lO^'' J m"^ 


Thermal electron momentum 


(p)e = e|Vp|« = 3.11T 


31 MeV 


Thermal photon momentum 


{p)vH = ell/pll = 2.1hT 


28 MeV 


Interparticle distance 


d^pl^-^l^ 


2.7 X 10-" m 


Coulomb Coupling Parameter 


A = eV(dT) 


5.3 X 10-=^ 


Effective Photon Mass 


rUyK = eT/3 


1 MeV 


Plasma Frequency 


uJpi — niph 


1.5 X 10^^ s"^ 


Dobye Screening Length 


Ad = l/iVSmph) 


1.1 X 10"" m 


EflFective Electron Mass 


rriF = er/(2v^) 

7e = e^T/{An) In(l/e) 

Fe = e*TV(37rs) ln(l/e) 

Fph = e*T'^/{GAT:E) \n{3.88E/e'^T) 


1.1 MeV 


Electron Damping Rate 


86 keV 


Electron Transport Rate 


0.54 keV for s = 19.3 


Photon Damping Rate 


0.70 keV for E = 2.75 T 


Electron Mean Free Path 




0.37 nrn 


Photon Mean Free Path 




0.28 nm 


Electron Collision Time 


Te = l/Fe 


1.2 X 10"^* s 


Photon Collision Time 


Tph ~ l/Fpfc 


9.4 X lO"^'' s 


Electron Viscosity 


Ve = 55.8TV[e^ln(l/e)] 
r^ph =48.7TV[e*ln(3.27/e)] 


7.9 X 10"' Pa s 


Photon Viscosity 


3.5 X 10^° Pa s 


Total Viscosity 


V = Ve+ Vph 


(1.1 - 1.6) X 10" Pas 


Electron Energy Loss 


dE/dx = e''T7(487r) ln(15.3S/e2T) 


200 MeV/nm for E = 100 
MeV 
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